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The recursion and path-integral methods are applied to analytically study the electronic structure 
of a neutral Ceo molecule. We employ a tight-binding Hamiltonian which considers both the s 
and p valence electrons of carbon. From the recursion method, we obtain closed-form analytic 
expressions for the 7r and a eigenvalues and eigenfunctions, including the highest occupied molecular 
orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) states, and the Green's 
functions. We also present the local densities of states around several ring clusters, which can be 
probed experimentally by using, for instance, a scanning tunneling microscope. From a path-integral 
method, identical results for the energy spectrum are also derived. In addition, the local density 
of states on one carbon atom is obtained; from this we can derive the degree of degeneracy of the 
energy levels. 

PACS numbers: 02.90.+p, 31.15.+q 



I. INTRODUCTION 

Since the recent discovery of a simple technique for the production in bulk quantities of fullerenes, undoped and 
doped C*6o molecules have generated enormous interest among chemists, physicists and materials scientists. In Cqq, 
carbon atoms sit at the sixty vertices of the pentagons and hexagons of a truncated-icosahcdron. The carbon-carbon 
bonds arc of two different lengths: 1.46 A for the single bonds (bonds on the pentagons) and 1.40 A for the double 
bonds (bonds on the hexagons not shared by a pentagon). The single (double) bonds are denoted by solid (dotted) 
lines throughout the figures in this paper. Ths-Jiterature on Cqq is vast, and here we do not attempt a review. The 
interested reader is referred to the referencesGrEI and papers listed therein. 

Most studies on the electronic structure of Cqo have been mainly numerical. In the present work, we focus on 
an analytical study of the electronic properties of a single neutral Cqq molecule. We model this system through a 
tight-binding Hamiltonian (Eq. (1)) which considers both the s and p valence electrons of carbon. This paper is 
organized as follows. In Sec. II, we discuss the physical nature and types of the interaction between these valence 
electrons, with a total number of 240. The couplings can be separated into two main contributions: coming from the 
7r-bonded and cr-bonded electrons. n 

We first focus on the it states. In Sec. Ill we use the recursion methodtl to analytically solve the Hamiltonian for 
the 7T electrons. We derive closed-form expressions for their energy eigenvaluess and eigenfunctions. The beauty of 
the recursion method for Cgo lies in the fact that the recurrence terminates very quickly {e.g., after four iterations 
only), providing exact and very concise expressions for the parameters of the recursion. The same results for the 
7r-states energy spectrum are derived in Sec. IV by using a path-integral method. This approach follows Feynman's 
programme: to compute physical quantities from sums over paths. In addition, the local density of states (LDOS) on 
one carbon atom is obtained; from this we can derive the degree of degeneracy of the energy levels. 

In the two sections mentioned above, every step of each approach is exact. All the calculations are done analytically, 
either by hand or with the assistance of computer symbolic manipulation software. Diagonalizations are fulfilled by 
iteratively applying the Hamiltonian on initial states. The relations between the recursion and path-integral methods 
are discussed in Sec. V. It can be seen that while these two approaches are related to each other, each one has its 
own advantages. These methods are significantly different from the ones used so far. Furthermore, they are neither 
numerical nor require the use of group theory. For a very lucid and clear account of a different approach, based on 
group theory and focused only on the 7r states, the reader is referred to Ref. [3]. 

In Sec. VI, we present the precise algebraic expressions as well as figures for the LDOS around several ring clusters 
(i.e., a carbon atom, a pentagon ring, and a hexagon ring; also for two opposite carbon atoms, pentagons, and 
hexagons). We find that around a pentagon (hexagon) ring the LDOS is large at low (high) energies. This is related 
to the fact that a pentagon (hexagon) has zero (three) double bonds and five (three) single bonds. These quantities 
are relevant to the several important experimental techniques which probe the local spectroscopy of molecules; (for 
instance, by using a scanning tunneling microscope, as described in the review in Ref. [9]). In Sec. VII we discuss the 
solutions of the c-states Hamiltonian Eq. (3). To describe the a states, the coupling between orbitals on the same 
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atom and that between, orbitals along the same bond are both taken into account. By following a "one band — two 
band" transformationt£l, the eigenvalues and eigenfunctions for the a states are then analytically obtained. 

II. ELECTRONIC STATES 

To investigate the electronic properties of a single Cqq molecule, we consider the four carbon valence electrons 2s, 
2p x , 2p y , and 2p z in the tight-binding Hamiltonian 

h t = 6a c «* Cm + E ta P c t a °iP- (!) 

i,cx <ij>,a,0 

Here, i denotes the carbon site and a denotes the valence orbitals 2s, 2p x , 2p y , and 2p z . Also, e a is the orbital 
energy and t a p is the hopping matrix element between orbitals on the nearest-neighboring sites i and j. The sixty 2p z 
orbitals, each pointing along the outward radial direction, are hybridized to form ir states. The other three orbitals 2s, 
2p x and 2p y , distributed on the plane tangential to the surface of the molecule, are hybridized along the lattice bonds 
to form a bonding and antibonding states. Let us assume that the mixture of these three orbitals at each site produce 
three sp 2 hybrid orbitals: sp 2 along the double bond and sp\ and sp 2 along the two single bonds, respectively. From 
a physical point of view, the 60 outer 7r orbitals are relevant to the conducting properties of the molecules and the 180 
a orbitals are mainly responsible for the elastic properties. The former are also responsible for the bond dimerization. 
Also, only 7r-states occur around the Fermi energy. 

As a result of its planar structure, the n electrons in graphite contain only pure 2p z orbitals. Nevertheless in the 
curved structure of Cqo, an extra small component from 2s, 2p x and 2p y is induced along the radial direction due to 
the surface curvature. The interaction between the 7r orbitals in Cqq are then slightly increased in comparison to that 
in graphite. However, since this component is very small, the overlap integral between n and a orbitals is negligible. 
The nonzero Hamiltonian parameters used here are the same ones as in Ref. [2]. Therefore, the original Hamiltonian 
(1) can be written in terms of H„. and Ho- as Ht = © H CT , where 

H* = -£ tijc+cj, (2) 

and 

i,a£P <ij> 

The electronic 7r-states Hamiltonian (2) describes the kinetic energy of the 2p z electrons hopping on the sixty vertices 
of a Cqo fullerene and Uj is the hopping integral between nearest-neighboring atoms i and j. The hopping integrals 
are considered unequal for single and double bonds. Also is the electron operator creating a 2p z orbital on the 
atom located at vertex i. A constant term corresponding to the on-site 2p z orbital energy is omitted in Eq. (2). In the 
cr-states Hamiltonian Eq. (3), a — sp 2 ,, sp\ and sp 2 denotes the hybridized orbitals. Also, V\ stands for the hopping 
integral between orbitals on the same carbon atom and V2 stands for that between orbitals on different atoms that 
are associated with the same bond (the length difference between single and double bonds is neglected here). A more 
detailed discussion on the origin of V\ and V2 is presented in Appendix A. 

Admittedly, this is a simplified model (like the "Ising model") for the electronic structure of Cqo, the "Hydrogen 
atom" of the main fullerene family of C 6 Q. n 2 molecules. However, its understanding is a convenient stepping stone 
to the study of more complex models. The spectroscopy of, for instance, the Hydrogen atom can be solved by 
using several different analytical approaches. These lead to the same analytical expressions for the eigenvalues and 
eigenvectors. In spite of the fact that modern computers can easily obtain numerical expressions for them, it is useful 
to have analytic results. In fact, analytical solutions have always been pursued in the area of spectroscopy of atoms, 
molecules, and clusters. In this paper, we pursue several approaches to the analytical study of the spectroscopy of an 
important molecule. We would like to emphasize that the approaches described here are not intended to substitute 
more traditional methods, but to present alternative viewpoints and complementary results. 

III. RECURSION METHOD APPROACH 
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A. Formulation 



We first focus on the tt states Hamiltonian (2). Let us begin with a brief outline of the recursion methodO for 
obtaining the eigenvalues and eigenfunctions. First, one must choose an appropriate normalized starting state |/o >• 
Further states are iteratively generated by the three-term recurrence relation 

H |/ n >= a n |/„ > +b n+1 |/„+i > +b n |/ n _i >, (4) 

with the condition &0I/-1 >= 0. Here, the real parameters a n and are determined by 

a n =< /„|H|/„ > (5) 

and 

b n+1 =< f n \n\f n+1 >=< /„+i|H|/„ > . (6) 

By convention, the &'s are chosen to be positive. The process terminates at |ujv_i > with the last recurrence 
H|/at-i >= aAr_ 1 |/ A r_i > +6jv-i|/jv-2 > • The constructed orthonormal states {|/o >j|/i >: -n - ,\In-i >} along 
with the parameters {ao, aj., . . . , <2jv-i} and {61, 62, ■ ■ ■ , bjv-i} thus constitute the "chain model"!! of a given Hamil- 
tonian. The representation of H in the new basis {|/„ >} is in a tridiagonal form. 

The energy levels can be achieved by constructing the following polynomials with P_i(E) = and Pq(E) = 1, 

p fj* (E - a n )Pn{E) — b n P n -i(E) 

P n +i(E) - . (7) 

Ora+l 

The eigenvalues E\ are determined by the N zeros of the last polynomial Pn(E) = 0, with an arbitrary nonzero bjy- 
It follows then that the eigenfunctions are 

1 rr ~ 1 

w Pn(Ex)\fn >, (8) 
A n=0 

where J\f\ — [J2n=o Pn(Ex)] 1 ^ 2 are the normalization constants for the eigenfunctions, and Pq(E\) always equals 1. 

So far, we have briefly outlined the recursion method. We now come to the application of this approach to the 
electronic structure of a fullerene molecule. For convenience, we work in units of the single-bond hopping-integral. 
The hopping amplitude is therefore 1 for every single bond and t for each double bond. It has been pointed outE that 
the value of t is about 1.1. We will use two alternative sets (denoted by A and 23) of starting states. Each set consists 
of two initial states from which the whole energy spectrum can be derived and the final results do not depend on the 
choice of initial bases. Both sets lead to the same solutions for the eigenvalues and eigenfunctions. From the physical 
point of view, this consistency ensures the equivalent interpretation of the results from these two sets of initial states. 
We denote by \j > the 2p z orbital centered at the Jth-atom. For convenience of visualization, in all our figures we 
flatten the truncated icosahedral structure of C^o into a plane. Note that labellings vary according to the different 
choices of initial states. 



B. Case A 



The first set of initial states starts from a 5-atom pentagon ring (see Fig. 1(a)) and a 6-atom hexagon ring (see 
Fig. 1(b)). Starting from a pentagon ring, we choose the initial state \uq > as a linear combination of the five orbitals 
on it 

1 5 

l M ° >= 7f E \i>- ( 9 ) 

Here the labelling is as shown in Fig. 1(a). From now on, H will donotc unless specified otherwise (e.g., H^). As 
an illustration, we present here the calculation of the first recurrence. 
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H |u > = [(|2 > +|5 > +t |6 >) + (|1 > +|3 > +t |7 >) 
V5 

+ (|2 > +|4 > +t |8 >) + (|3 > +|5 > +t |9 >) + (|1 > +|4 > +t 1 10 >)] 
= -2 [-L (|i > +|2 > +|3 > +|4 > +|5 >)] +* (|6 > +|7 > +|8 > +|9 > +|10 >)] 

= -2 \u Q > +t \ui > . 

We then obtain a = —2, b\ = t and |ui >= — J2]=6 \j > • Following the same procedure, we can construct further 
states and obtain the parameters a's and 6's at the same time. The beauty of the recursion method for 6*60 lies in 
the fact that the recurrence terminates very quickly, exactly at \ur >. As a result, we have exact and very concise 
formulas for ao through 07, and 61 through 67. Similarly, starting from a linear combination of the six orbitals on a 
hexagon ring, we have the initial state \vo > as 

\vo>=±p^-l) i+1 \j>, (10) 

where the labelling is referred to Fig. 1(b). The choice of alternating signs, (— meets the three-fold symmetry 
requirement of a hexagon ring in Cqq, with alternating single and double bonds. Again, the recurrence terminates at 
\vj >. We then obtain another group of parameters a's and 6's. All these \u n > and \v n > states, as well as their 
respective parameters a n and b n +i, are summarized in Tabic I. 

Through these two groups of a's and 6's, we can respectively construct two polynomials P£(E) and P${E) which 
are of eighth degree in E. Here the superscript p, for pentagon, (h for hexagon) refers to the polynomial constructed 
from the group of a's and 6's generated by |u > (l w o >)• The roots of these two polynomials can be analytically 
obtained and correspond to the electronic energy levels. It follows then that we have a total of 16 distinct eigenvalues. 
When t = 1, a common root —1 exists for both polynomials Pg and P^. We thus have 15 distinct energy levels. It is 
also straightforward to obtain the eigenvectors through Eq. (6). 



C. Case B 



In this approach, we exploit the symmetry property that the inversion operator leaves the C@o molecule invariant. 
We therefore take the first (second) starting state \(j>o > (\ipo >) as a linear combination of the orbitals on two opposite 
(i.e., antipodes) pentagon (hexagon) rings 

and 

1^ >= ~m i { - iy+1 [|j > +v[i ' >]> (12) 

where the site labels are presented in Fig. 1(a) (1(b)) for the first (second) starting state. Here V stands for parity, 
with the value +1 or —1. Note that atom j and atom f are antipodes. Following the same procedure, we find the very 
convenient and remarkable result that the recurrence terminates even faster at \<p3 > Qtpg >) for the first (second) 
initial state. Due to the two possible values of V, we have two different results for a 3 for each starting state. In 
Appendix B we present the calculation of the last recurrence, namely those for \<j) 3 > and |<^ 3 >. It can be clearly 
seen that H \ (j> 3 > generates two possible values for 03 and no further state; and so does H | ^3 >. As a consequence, 
we obtain four groups of {00,01,02,03} and {61,62,63}. All the states \<p n > and \<p n >, as well as their respective 
parameters a n and 6„ + i are listed in Table II. It follows that we can construct four polynomials Pl + (E), P%~(E), 
P£ + (E) and P£~(E). Each one is fourth degree in E. The superscript p+ stands for the polynomial constructed 
from |(/) > with V equal to +1, and similarly for the others. With the choice of 64 = I/6162&3, these polynomials can 
be explicitly written as 



Pf + {E) = E 4 + (2t + 3) E 3 + (5t - 1) E 2 - (2t 3 - t 2 + 8) E - [t + 2){t 3 -t 2 +t + 2), 
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Pl~{E) = E 4 + 3E 3 - (2t 2 - t + 1) E 2 - (3t 2 - At + 8) E + (t A - t 3 + t 2 + At - 4), 



P%+(E) = E 4 - 2 (t + 1) E 3 + (3 1 - 1) E 2 - (2 1 3 + t + 2) E - (t 2 + 1) (t 2 + t - 1), 

and 

P^(E) = E i - 2E 3 - (2t 2 + t + 1)E 2 - (2t 2 +t + 2)E-(t + l) 2 (t 2 - t + 1). 

By analytically solving these four polynomials, we obtain the same 16 eigenvalues obtained above (case A). Similarly 
we can obtain the eigenvectors, which are also equal to those obtained from the alternative A. 



D. Results 



In Table III, we summarize the eigenvalues and the corresponding eigenvectors for the case t = 1. In Table IV, we 
present the closed-form eigenvalues explicitly expressed in terms of the single-bond hopping integral t\ and the double- 
bond hopping integral t-2.. Thus, eigenvalues for the limiting cases t\ = or t2 — can also be readily inferred. As to 
the cigenfunctions, we present those for the highest occupied molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO) in terms of t, the ratio of the hopping matrix elements for single and double bonds. Since 
the respective degree of degeneracy for each eigenvalue can be acquired from the local density of states on a carbon 
atom (discussed in a later section), the other degenerate eigenvectors can be generated by standard group theory 
analysis. This is outside the scope of the present paper. We therefore only present these eigenfunctions derived from 
the pure application of the recursion method. It is instructive to notice that the following points further support the 
results for the degeneracy. They are: (1) the sum of the product of energy and its corresponding egeneracy equals 
the trace of the Hamiltonian, which is zero; (2) the number of states with even parity equals that with odd parity, 
and (3) the behavior of the eigenvalues can be easily studied in the limits when either t\ = or t 2 = 0. 



E. Relations between alternatives A and B 



It is worthwhile to point out the following relations between the cases A and B presented before. First, Pl{E) = 
Pl + {E)Pl~(E) and Pg{E) = P% + (E) P^ (E) up to an overall constant factor. In other words, roots solved from 
Pl(E) (Pg{E)) are identical to those solved from P% + (E) and P%~ (E) (P% + {E) and P£~(E)). Second, from Tables 
I and II we can see that, for n — 0, 1, 2 and 3, 



: (\u n > -V\u 7 - n 



\4>n >= — (\U„ > -V\u 7 -n >) 

and 

\<Pn >= -j= (K > +V\v7-n >)■ 

Third, for those eigenvalues E\ which are common roots of PP and P 4 P+ , or common roots of P£ and P£~ , 

P 7 - n (E X ) = ~P n (E X ). 

Also for those eigenvalues E\ which are common roots of Pg and P][~ , or common roots of P$ and P£ + , 

P 7 - n (E x ) = P n (E x ). 

Here, all the P n (E\)'s refer to the polynomials constructed in the case A and n = 0,1,2 and 3. Fourth, Pi(E\), 
P2(E\) and P3(E\) calculated from the two alternatives A and B are the same. Fifth, for the eigenvector with 
respect to the same eigenvalue, the normalization constants AOi and Mb, calculated in A and B respectively, satisfy 
Na = \/2 A/g. From the above five properties, the equivalence of results from both alternatives becomes clear. As an 
example, for an eigenvalue E\ which is a common root of and P% + , we can obtain the eigenvector >from A 
as 
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I** > = AT E P n(Ex)\Un > 
n— 

= ac* [P ° 71 > ~' M7 >) + Pl 71 > ~'" 6 >) 

+P2 71 (|U2 > _|us >) + Pa 71 (l " 3 > _l " 4 >)] 
i 3 

= W E P n(^A)|^„ >, 
^ n— 

where the last equality is just the result obtained directly from B. In fact, from the results of alternative A, we can 
understand the parity property associated with Cqq. 



IV. PATH-INTEGRAL (OR MOMENT) APPROACH 
A. Methodology and application 

In this method, the first and central task is the computation of path-integrals (moments), defined by 

M l =<m i \j>, (13) 

where the order I is a positive integer. Note that Mo = 1 when I = 0. The physical meaning of the above quantum- 
mechanical expectation value is as follows. The Hamiltonian H is applied I times to an initial 2p z electron state \j >, 
localized at carbon site j. Each time H is applied, the electron gains a certain amount of kinetic energy depending 
upon the bond (single or double) it travels. This enables the electron to hop through I bonds, reaching the final 
state H'|j >. The path- integral (moment) just equals the total kinetic energy gained by the electron returning to 
the starting site j after hopping I steps. It is obvious that Mi will be zero when the Z-hops path does not return to 
the starting site. In other words, Mi = when there is no path of I hops for which the electron may return to the 
initial site. For the case t = 1, the absolute value of the moment Mi is the total number of closed paths of Z-steps 
starting and ending at the same site. This approach follows Feynman's programme: to compute physical quantities 
from sums over paths. 

The path-integrals can be calculated analytically by hand, as well as by computer using a symbolic manipulation 
program. Below we describe both ways. Let us first examine the action of the Hamiltonian on an arbitrary state 
\j >. This results in three nearest-neighbor atom states with an additional factor accounting for the respective bond 
hopping energy. For example, H |1 >= (— t) |2 > +(—1) |3 > +(—1) |3" > (labels as shown in Fig. 2). For simplicity 
and without any loss of generality, we choose orbital |1 > to be our starting state. Now, for 1 = 1, starting from 
vertex 1 and following the connectivity of Cqq (Fig. 2), we write down the factors —t, —1 and —1 on the vertices 2, 
3 and 3" respectively. A similar procedure holds for 1 = 2. Starting from the three resulting vertices with respective 
factors, we then write down (— 1) = t on the vertices 4, 4', 4" and 4"'; (— 1)(— 1) = 1 on the vertices 5 and 5"; and 
(— 1)(— 1) + (— 1)(— 1) + (— t) = 2 + i 2 on the vertex 1. Our strategy is: each time the power of the Hamiltonian (i.e. 
kinetic energy) increases by one, we move to the adjacent nearest-neighboring vertices. Also the factor on each vertex 
is the sum of the factors on the three nearest-neighbor vertices times the bond hopping integral between vertices. It 
is then straightforward, with the aid of "the flattened 6*60 graph" , to obtain all the factors on the available vertices 
for any power of the Hamiltonian. From its definition, it is evident that the path-integral of order I is just the factor 
on the vertex 1 for H*. For example, Mi = (i.e., the number of closed paths obtained by moving one step is zero) 
and M2 = 2 + i 2 . By following this strategy, we can generate, one by one, all the path-integrals to any order. It is 
worthwhile to notice the mirror-symmetry between the left and right halves of Cqo for the atom vertices and bonds. 
So we need only concern ourselves with the factors on the vertices in the right part. Furthermore, because of the 
geometrical equivalence, with respect to vertex 1, of eight pairs of vertices (j and f in Fig. 2), the total number of 
independent vertices can be reduced to 24. All calculations up to this point can be done analytically by hand. 

An alternative analogous procedure pan be implemented by using symbolic manipulation software on a computer. 
First, we define an auxiliary quantityU, Wi(j), which is the sum over all possible paths of I steps on which an 
electron may hop from the vertex 1 to the vertex j. From the connectivity of Ceo, we can then construct 24 
independent recurrence relations. For instance, Wj + i(l) = -tWj(2) - 2Wj(3), Wj+i(2) = -tWi(l) - 2Wj(4), and 
Wa+i(3) = —tWi(A) — Wi(l) — Wi(h). The recurrence relations state that the vertex j can be reached by taking the 
(Z+l)th step from the three nearest-neighbor vertices. The factors —t and —1 account for the connecting bond-hopping 
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integral. With the initial conditions Wrj(l) — 1 and Wo(j) — for the rest of j's, we can obtain the path-integrals to 
any order as Mi = Wi(l). We list M 3 through M 8 here: 0, t 4 + 8t 2 + 6 (15), -2, t 6 + 18 i 4 + 4t 3 + 48 t 2 + 20 
(91), -14 1 2 - 14 (-28), t 8 + 32i 6 + 16t 5 + 184/; 4 + 48i 3 + 256i 2 + 70 (607), for I = 3, ■ • • , 8. The numbers in 
parentheses are the values when t = 1. The correctness of the calculated path- integrals is assured by the consistency 
of the results from these two approaches. It is evident that through these two approaches we can also obtain the 
quantities < 1|H \j > for j ^ 1 which can be appropriately interpreted as the "sum-over-paths'" between sites 1 and 
j. For instance, < l|H'|j > just equals Wi(j). 

To obtain the energy spectrum, we again utilize equation (7). So, the next step is to express the parameters a n and 
b n +i in terms of the moments. We employ the following formulas^: define the auxiliary matrix M with the first row 
elements defined as Mgj = Mi. The other rows are evaluated by using only one immediate predecessor row: 

M n .i = 77 772 / ^ M n,k M n -x,i-k , n> l;i = 0, I,---. (14) 

Mn-1,2 - M„_ 11 ^ 

The a„'s and 6 n +i's are obtained from the elements of the second and third columns as 

a„ = M n ,i (15) 

and 

bn+i = M n< 2 — M 2 X , (16) 

here n = 0, 1, 2, • ■ ■ . Note that elements in the first column M n> o are always equal to 1. We analytically find that 
b 2 5 exactly equals for t — 1 and 6f 6 exactly equals for an arbitrary t. Below we discuss the t = 1 case. Similar 
conclusions can be drawn for t ^ 1. The appearance of absolute zero for 6f 5 in the case t = 1 indicates the truncation 
at 615. In this sense, we expect 15 eigenvalues. It also turns out that the highest order of moment we need is 
I = 30. The moments Mg through M30 for t = 1 are presented in Table VI. Through the calculated parameters 
{do, 01, . . . , 014} and {61, 62, ... , 614} (listed in Table VII), we can construct the polynomial P\§(E) by using Eq. (5). 
By solving P\§{E) = 0, we thus obtain 15 energy levels. The results are exactly identical to those obtained from the 
recursion method approaches. We present Mg through Mw for an arbitrary t in Table VIII. 



B. Alternative application of the path-integral method 

In the above description of the path-integral approach, the main ingredient is the computation of the path-integrals 
for a natural choice of state |1 > centered at the atom labelled by 1. However, it is worthwhile to incorporate the 
inversion symmetry property. Therefore, instead of focusing on a single localized state, we turn to the computation 
of the path-integrals with respect to states 

\I± >= -j= (|1 > ±|24 >), (17) 

where atoms 1 and 24 are antipodes (as labelled in Fig. 2). It is a simple exercise to construct the following identity 
for the path-integrals defined by < l±|H'|/± >, 

< /±|H'|/± >=< 1|H'|1 > ± < 1|H'|24 > . (18) 

These moments can then be readily obtained since the quantities < 1|H'|1 > and < 1|H'|24 > are already available. 
Note that the lowest order for the appearance of a nonzero < 1|H Z |24 > is I — 9. This is because the shortest path 
for vertex 1 to reach vertex 24 contains 9 steps. We thus obviously have < /±|H'|7± >= Mi for I < 8. 

For the case t = 1, the highest order we need for < 2 + |H'|i + > (< 2_|H J |J_ >) is I = 14 (16), because b 2 
(6g) calculated from moments < i+|H'|i+ > (< J_|H |i_ >) turns out to be exactly 0. Through the computed 
parameters {arj, . . . , a§} and {61, . . . , b$} ({ag, ■ ■ ■ , &7} and {b\, . . . , 67}), we analytically obtain 7 (8) eigenvalues which 
are identical to those belonging to the V = +1 (V = —1) category from the recursion method. Moments < i+|H i |i + > 
(I = 9, • • • , 14) and < 7_|H |i_ > (/ = 9, • • • , 16) are presented in Table IX and their respective set of parameters 
a n 's and 6„+i's in Table X. The result that a n and 6 n +i for n — 0, 1,2 and 3 in Table X are identical to those in 
Table VII comes from the fact < 7±|H'|/± >= Mi for I < 8. For an arbitrary t, we need moments < 7±|H'|/± > up 
to order 16 ( < I±|H 9 |i± > through < 7±|H 16 |/ ± > are listed in Table XI). We analytically find that 6| calculated 
from these two sets of moments exactly equals in both cases. Consequently, consistent results for the eigenvalues 
are recovered. In this section, we have presented an unconventional choice of initial states and concentrated on the 
moments with respect to these states. It is shown that this approach is even more efficient in analytically obtaining 
the energy eigenvalues. 
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V. RELATIONSHIP BETWEEN THE RECURSION AND PATH-INTEGRAL METHODS 



Generally speaking, the path-integral method is closely related to the recursion method, especially in the aspect 
that both methods lead to the same results for the parameters a n and b^ +1 . In this section, we illuminate this point 
by showing that the same expressions for the parameters in the recursion method case B can be obtained through the 
moment method. As the initial states are \<pQ > and \ipo > in this case, the moments we now need to compute are 
< <^o|H'|0o > and < </?o|H'|</?o >• It is straightforward to find that 

< c^olH^o > = < 1|H'|1 > +2 < 1|H'|3 > +2 < 1|H'|5 > 

+ V{< 1|H'|24 > +2 < 1|H'|22 > +2 < 1|H ; |19 >), (19) 

and 

< (p Q \H l \ip > = < l\H l \l > - < l\H l \2 > - < l\H l \3 > - < l\H l \7 > +2 < 1|H'|4 > + 
+ V (< 1|H'|24 > - < 1|H'|23 > - < 1|H'|22 > 

- < 1|H'|21 > +2 < 1|H'|20 >). (20) 

Notice that the site labels on the right hand sides of the above two equations refer to Fig. 2. Through the techniques 
for the calculation of these quantities previously discussed, the moments < (fo|H'|<^o > and < (^o|H'|<^o > can be 
readily obtained. 

Anticipating the termination at b 2 , we only need these moments up to order 8. In Table XII we give the moments 
with respect to |</>o > and |<^o >. By utilizing Eqs. (14-16), we thus obtain a„'s and 6^ + i's which are consistent with 
those (in Table II) derived directly from the recursion method. To illustrate this consistency, we explicitly present 
here the results for the b£ +1 's. From < O |H'|0 O >, we have b\ = t 2 , bl = 2, b 2 = 1, and b\ = (1 - V 2 ) t 2 = 0. 
From < <y9o|H'|(^o >, we have b\ = 1, b\ = t 2 , b 2 = 1, and b\ = (1 — V 2 ) t 2 = 0. We thus demonstrate the fact that 
the same results for the parameters a n and £>^ +1 can be obtained by using either the recursion or moment methods. 
The advantage of the recursion method lies in that we can simultaneously generate the states and the parameters. 
However, it is sometimes difficult to derive the states and parameters when the recursion method is applied to some 
starting state, for example, a single carbon atom state \j >, while the moment method provides standard procedures 
to calculate the parameters after the moments are obtained. 



VI. LOCAL DENSITY OF STATES 



In the recursion a»d-momcnt methods, the diagonal element of the Green function (E — H) 1 can be expressed as 
a continued fractiontrtHl 



Go(£) = </o|-^jj|/o> 



b{ 

E-a Q -J 



(21) 



E — a\ 



E - 02 - • 



E — ajv-i 

The local density of states p(E) for |/o > is related to the imaginary part of Gq(E) by 

p(E) = lim — Im G (E + ie). (22) 

From the computational point of view, Go (E) can be obtained by iteratively applying the following transformation 

1 



G n {E) - ■ - 7 



E — a n — b n+1 G n +i(E) ' 



starting from Gn^\{E) — 1/ (E— ajv-i). By plugging the parameters a n and b\ +1 from either one of the two methods 
into Eq. (21) and using Eq. (22), we thus obtain the local densities of states on several initial states. In our two 
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alternative applications of the recursion method, we have used four different starting states: \uq >, \vq >, |0o > and 
\ipo >■ Their Gq(E)'s can be written explicitly as 

G^ 1 = E + 2 - t 2 {E - 2 [E + t - (E + 1 + GpY 1 ]' 1 }' 1 (23) 

and 

Gq 1 = E - t - 1 - {E - t 2 [E - 1 - (E - Gh)- 1 }- 1 }' 1 , (24) 

where 

r _ J V t for the state 1 4>o > 

= \ -t 2 {E + 1 - [E + t - 2 (E - t 2 (E + 2)- 1 )- 1 ]- 1 }- 1 for the state \u > 

and 

c { Vt for the state |<^o > 

^ h ~ \ t 2 {E - [E-l-t 2 (E- (E-t- lr 1 )" 1 ]- 1 }- 1 for the state \v >■ 

The LDOS around these ring-clusters are plotted in Fig. 3 and Fig. 4. Notice that around a pentagon (hexagon) 
ring the LDOS is large at low (high) energies. This is related to the fact that a pentagon (hexagon) has zero (three) 
double bonds and five (three) single bonds. Also from the path-integral approach, the local densities of states on 
two antipode carbon atoms are plotted in Fig. 5. In principle, they are experimentally accessible by using a scanning 
tunneling microscopeB. In one of the path-integral method approaches, the initial state is a 2p z orbital on a carbon 
atom. The local density of states in the case t = 1 is plotted in Fig. 6. From Fig. 6, we can obtain the degree of 
degeneracy for each energy level which is the respective LDOS value times 2. The common factor of 2 comes from the 
consideration that the total number of tt electrons is 60. For t ^ 1, conclusions for the degeneracy can be similarly 
drawn. 



VII. SOLUTION FOR THE ELECTRONIC a STATES 



It is apparent that the Hamiltonian H a is more complex than H ff due to the two different couplings. However, 
an analytic transformation of this. Hamiltonian into a simpler one with a single "renormalized" hopping parameter 
between sites can be establishedtL3. The energy states can then be readily solved. The details of this transformation 
as well as the subsequent solutions for the a states are given in Appendix C. As a result, we obtain 90 bonding a 
states and 90 antibonding a states. Among the bonding (antibonding) states, 30 states are lumped together at the 
energy level V1 — V2 (Vi + V2). The other 60 a bonding and 60 a antibonding states are closely related to the energy 
spectrum for the 7r states. To be more specific, the eigenvalues of the other 60 bonding (designated by plus sign) and 
60 antibonding a states (designated by minus sign) are 

" ^ (25) 

Avr (25) 

where E\ are the -k state eigenvalues listed in Table III. 




VIII. CONCLUSION 



In conclusion, we use several approaches based on the recursion and path-integral methods in order to explore the 
electronic structure of a Cqq molecule, obtaining exact closed-form expressions for the for the tt and a eigenvalues 
and eigenfunctions, including the HOMO and LUMO states, as well as the Green's functions and LDOS through 
alternative methods. These quantities are relevant to the several important experimental techniques which probe the 
local spectroscopy of molecules (for instance, by using a scanning tunneling microscope, as described in the review 
in Ref. [8]). For comparison purposes, we have also done a direct numerical diagonalization of the full Hamiltonian 
and the results obtained are consistent with those from the previous analytical methods. However, the much more 
elegant and powerful recursion and moment methods provide valuable insights and closed-form expressions for various 
aspects of the electronic structure of Geo • 
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APPENDIX A: ORIGIN OF Vi AND V 2 



In this appendix we examine more closely the physical interpretation for the hopping integrals V\ and V 2 . Let 
i, a > be a hybrid orbital located at site i. We have 

\i, sp 2 a >= a\ \2s > +4 \2p x > +4 \2p y >, (Al) 



\i, sp 2 b >= b\ \2s > +b l x \2p x > +b\ \2 Py >, (A2) 

and 

\j, spl >= a{ \2s > +< \2p x > +a{ \2 Py >, (A3) 

where i and j are nearest neighbors. We also assume that orbitals \i, spl > and \j, spl > lie along the bond connecting 
i and j. It is straightforward to obtain 

- Vi =< i, spl\K a \i, spl >= a\*b\ e s + (4% + a L y *b y ) e p , (A4) 

and 

- V 2 = < i,sp 2 a \H a \j,spl > 

= 4*4 t ss + a^ai t PxPx + a y *o? y t PyPy 

+ (4*4 + 4*4) t sPx + (4*4 + 4*4) t sPy + (4*4 + 4*4) t PxPy . (A5) 

Here e s and e p are the 2s-level and 2p-level energies and the i's are hopping matrix elements between orbitals on 
nearest-neighboring sites. It is evident that, with a suitable choice of a local coordinate system, we can always have 

< i, sp 2 a \H a \i, spl >=< i,sp 2 a \Tl a \i,sp 2 c >=< i,spl\R a \i,sp 2 c > (A6) 

and 

< hsp 2 a \U a \j,spl >=< i, spl\H a \k, spl >=< i,spl\U a \l,spl > . (A7) 

In the above equation, we have assumed that orbitals \i, spl >, \k, spl > lie along the bond connecting i and k and 
orbitals \i, spl >, \l, spl > along the bond connecting i and I. 



APPENDIX B: CALCULATION OF H|0 3 > AND H \ip 3 > 



In this appendix we present the calculation of the recurrence relations for states |</>3 > and \ip3 > in the recursion 
method case B. It is shown that the recurrence terminates at these states and two solutions for (%■$ are constructed 
from each state. 

^20 ^30 30 



H i& > = -7= E (Ij > w >) + ~m [ E (|j > +t b " >} + r E {t |j > +l/ >)] 

j=ll .7=21 i=21 

1 30 

>-{--= E [(i + ^*)li>+(* + ^)l/>]} 



20 ,=21 



30 



\4> 2> -{-^= <r [(i+vt)\j>+(i+vt)v\j>>]} 

V 20 ,=21 

|^2 > -(l+7>t) |^3 > • (Bl) 
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H|y»3 > = |v>2 > +-^=[t(|23 / > -|24' > +|26' > -|27' > +|29' > -|30' >) 



+ Vt{\23 > -|24 > +|26 > -|27 > +|29 > — 130 >)] 

\<P2 > +Vt-^=[ (|23 > -|24 > +|26 > -|27 > +|29 > — j30 >) 
v 12 



+ "P(|23' > -| 24' > +|26' > -| 27' > +|29' > -|30' >)] 
= |(^ 2 > +Vt \tp 3 > . (B2) 



APPENDIX C: ONE BAND-TWO BAND TRANSFORMATION 

In this appendix we derive the transformation of the tr-statcs Hamiltonian into a 7r-type one. The resulting solutions 
of the eigenvalues for the a states are also presented. We first write the a states eigenfunctions as 

i 

where i runs over the atoms and C = a,b and c represents the corresponding coefficient for the orbital a at each 
site. Let us now focus attention on a given atom labelled i. Also let atoms j, k and I be the three adjacent nearest 
neighbors of i-atom. The energy eigenvalue equation H CT 1^ >= e |\I> CT > then reads: 

ea i = -V 2 a j -V 1 {h + a), (C2) 
e bi = -V 2 c fe - Vi {a, + a), (C3) 



and 



e = -V 2 b t - V x (a t + bi), (C4) 

e cij = -V 2 ftj - Vi (bj +Cj), (C5) 
e c k = -V 2 bi - V x (a k + b k ), (C6) 



e 



= -^2 - 1^ (a t + c,), (C7) 



and so on. We then define 



Ai=ai + bi + Ci, (C8) 

as the sum of the coefficients for orbitals at the site i, and 

Bi = aj+c k + bi, (C9) 

as the sum of the coefficients for orbitals on the bonds associated with the site i. Similar equations hold for the other 
A's and B's. In the first place, Eqs. (C2) and (C5) can be rewritten as 

{e-V 1 )a i + V 2 a j = -V 1 A i , (CIO) 

and 

V 2 <n + (e-Vi) a, =-Vi Aj. (Cll) 
When the situation Aj, = Aj = occurs, the secular equation immediately yields 

e = V 1 ±V 2 . (C12) 
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In all, there are thirty pairs of similar equation sets. We thus have two energy levels V\ — V 2 belonging to the bonding 
states and V\ + V 2 belonging to the antibonding states, each having a degree of degeneracy equal to 30. Secondly, by 
respectively summing up Eqs. (C2) through (C4) and Eqs. (C5) through (C7), we have 

eA i = -V 2 B i -2V 1 A i , (C13) 

and 

(e - Vi) B t = -V 2 A, - V Y (Aj + A k + A t ). (C14) 

By substituting Eq. (C13) into Eq. (C14), we obtain 

[(e + 2 Vi)(e - V{) - Vi] A l = Vi V 2 (Aj + A k + At). (C15) 

It can be directly recognized that Eq. (C15) is entirely equivalent to the problem of a tight-binding Hamiltonian with 
one state per atom and a single nearest-neighbor hopping integral, whose eigenvalue and eigenvector solutions are 
already fully explored in the text. Let E\ (A = 1, 2, • • • , 15) stand for the eigenvalues summarized in Table III. It 
follows then that 

(e x + 2 V^ex - Vi) - V? = Vi V 2 E x . (C16) 

Thus, we obtain 



ex = ~ ± \ + 2 Vi + V 2 2 + V 1 V 2 Ex 




1 + v 2 Ex + IW- (C17) 



The minus sign designates the bonding states, and the plus sign designates the antibonding states. Finally, in the 
limit V 2 » Vi, we have 

V Vi 

ex - - y ± V 2 ± -1 Ex. (C18) 

It is also straightward to obtain the corresponding eigenvector for an eigenvalue ex- Since all the Ai's are already 
known, the initial coefficients ai,bi and Cj can be calculated. 
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TABLE I. States and parameters for case A in the recursion method approach. It is interesting to notice that a n — a 7 _ 
and b n — bs- n for n = 4, 5, 6, 7. 





Starting from a pentagon ring 


Starting from a hexagon ring 






n 


\u„ > a n b n 

1 • . ^ 


\Vn > 

1 b / h \ i -4- 1 l - . 


a„ 


b n 



b EL \J > 



4=Y™ \i> 

/TO ^3 = 11 U 

E?=2i b > 



1 E-° 



b"> 

n 1/ > 



vTo 

T7f Ej= 6 b'' > 

e;=i b' > 





-1 
-1 
-i 



-2 



1 
1 

t 



V6 



^15 



4se: 7 (-i) j b> 



x/6 E. 



, , 4 (-1) J+1 [|J >+b' + 3>+|i + 6>] 

24 / t\j + lr 



(-1) J+ [|j >+|j + 3>+|j + 6>] 



7S E - 4 =23 (-l) J+1 [b" > +10' + 3)' > + 6)' >] 
75 E-Li4 (-l) i+1 [b-' > +10' + 3)' > +10 + 6)' >] 

^ E" 7 (- W > 



s/6 



EL (-iF +1 b" > 



o 
i 




1 



1 + t 



TABLE II. States and parameters for case B in the recursion method approach. The parity V can be +1 or — 1. 





Starting from two opposite pentagon rings 






n 


\4>n > 


a n 


b n 





t% e; =1 lb > +r \f >] 


-2 




1 


-tttjE" [b>+^b" >] 





t 


2 


TfoE-L [\J>+V\f >] 


-t 




3 




-(1+Vt) 


1 




Starting from two opposite hexagon rings 






n 


\lfin > 


a„ 


bn 





tbEh (-i) 3+1 [b>+^b" >] 


l + t 




1 


715 E" (-l) J [|i>+P b" >] 





1 


2 


7% Ej-14 > +IJ + 3 > + |j + 6> +7>(|j' > +10 + 3)' > +10 + 6)' >)] 


1 


t 


3 


7T5 E - 4 =23 (-l)' +1 [b > +b' + 3 > +\J + 6 > +P (\f > +10 + 3)' > +10 + 6)' >)] 


Vt 


1 
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TABLE III. The eigenvalues E\ and the corresponding eigenvectors X)n=o Pn(E\)\f„ > for Ceo with t = 1. Here > 
j((9* >) denotes \<j> n > (\<p n >) with V = ±1. 



P2(E X ) 

-V2 
^2(3 + x/l3)/4 
v/2(3- VT3)/4 



V , 2[6+2y / 5+(y / 5+l) a] 
16 

V^[6+2y / 5-(v / 5+l) q] 



-V2 
V2/2 
y/2/2 



V2[-4+8v^+(v / 5-l) a] 
16 

v^[-4+8v / 5-(y / 5-l) a] 
16 

y^[-4-8v^-( N /5+l)/;] 
16 

V2l-4-8VT, + (V5+l)l3] 
16 



-1 
-1 

(3 - v / 5)/2 
(3 + V5)/2 

(1 - %/5)/2 
(l + V5)/2 

-4 + VTf 

-4-y/Vf 



l/n > 



-3 
-1 

(-1 + x/l3)/2 
(-1- Vl3)/2 

(-3 + V5 + a)/4 
(-3 + x/5-a)/4 
(-3- V5 + p)/4 
-3- V5-/3)/4 



(3 + >/5)/2 
(3 - x/5)/2 

(1 + V5)/2 
(1 - a/5)/2 
(l + VT7)/2 
(1-x/17)/2 



(3 + x/T3)/2 
(3 - x/T3)/2 



5+%/5+a 
4 

5-j- y_5— q 

4 

5-x/5+/3 
4 

5-V5-/3 





-3 

(-1 + v/5)/2 
(-1 - x/5/2 

(-3 + v / 5)/2 
(-3- V5)/2 
(-3 + vT7)/2 
(-3 - vT7)/2 



16 

v^[6-2v / 5-(y/5-l) g] 
16 

y / 2[6-2Vli+(v / 5-l) ) 3] 
16 



-1 

2 

(-1 + 71)12 
(-1 - 71)12 

(-1 - x/5)/2 
(-1 + y/§)/2 
(5 - vTf)/2 
(5 + vTT)/2 



(39 + 9\/T3)/4 
(39 - 9x/l3)/4 

95+5^ _|_ (25+75) a 



16 



95+5^5 _ (25+\/5) a 
8 16 

95-5x/5 , (25-y/5)|9 
8 16 

95-5^5 _ (25-VS)J 
8 16 



3 
15 

(15-5\/5)/2 
(15 + 5\/5)/2 

(15-3V5)/2 
(15 + 3VE)/2 
51 - 12^17 
51 + 12^17 



Here a = y/ 2(19 + 75) and (3 = y/ 2(19 - 75). 
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TABLE IV. The eigenvalues and the corresponding degree 
of degeneracy for Ceo with a single-bond hopping-integral ti 
and a double-bond hopping-integral ti. The characteristic 
polynomials from which those eigenvalues are solved are indi- 
cated in the left-hand column. 



pv 



pP- 
M 



Energy 

-(2*i+< 2 ) 
-(ti+i 2 )/3 + 2f7 
-(*i + * 2 ) /3 - 77 + a/3£ 
-(ti+ti)/3-r)-V3(, 

[(-3 + VE)t 1 +r] /4 
[(-3 + >/5)ti-T] /4 
[(-3- \/5)*i+ 7 ] /4 
[(-3-a/5)*i- 7 ] /4 

(*i + a /5* 2 + 4* 2 ) /2 

(ti - yoTf+Ttf) /2 

* 2 + (1 + a/5)*i/2 
* 2 + (1-a/5)*i/2 



Degeneracy 
1 
5 
5 
5 

3 
3 
3 
3 

4 
4 
3 
3 



(ti + a /5* 2 -4*! * 2 +4*| ) /2 5 

ph- (*i - y /5tf-4tit 2 +4t| ) /2 5 

4 (*i + y /5* 2 + 8*i* 2 +4* 2 ) /2 4 

(*i - v/5t?+8ti * 2 +4*|) /2 4 

Here r = a/2 [8*| - (4 + 4a/5)*2*i + (15 + 5>/5) *?] and 7 = 
a/2 [ 8 t\ - (4 - 4 a/5) ti *i + (15 - 5^5) *fj. Also jj and £ sat- 
is fy (16*1 ~ 24ii*i + 12_Mi + 25*?) /54 = ?? 3 - 3^ 2 and 
a/3 (64 *| - 160 *| *i + 288 *| *f - 200 * 2 *f + 125 tf) /18 
3^-C 3 - 



TABLE V. The HOMO wave function ^- £^ =0 p »( s a)|^ > and LUMO wave function ^ 
in terms of *. Here \ip„ > stands for \tp n > with V = — 1 and |0~ > for \4> n > with "P = — 1. 



=0 Pn(i?A 



> for Ceo 



HOMO 



Pi(E x ) 
Pi(E x ) 
Ps(Ex) 
A/a 



(1 - a/4* 2 -4* + 5)/2 



6* 2 



-*- (I + a/4* 2 
*- (3 
*-(l 
6* 



4* + 5 )/2 
a /4* 2 -4*+~5 )/2 
a/4* 2 -4* + 5)/ 2 



f + (3* 



)a/4* 2 -4* + 5 



LUMO 



E x 
Pi(Ex) 

P2{Ex) 

P 3 (Ex) 
Nx 



(-3 -a/5 + 7) /4 
_ (5-a/5 + 7 )/4* 
a/2 [ 10 - 6 a/5 - 4 (1 - a/5) * + (1 - a/5) 7]/ 16 * 
-a/2 [4 a/5 + 4 (1 + a/5) * + (1 + a/5) 7 ]/ 16 * 
[ 150 - 50 a/5 + 80 * 2 - 40 (1 - a/5) * + (25 - 5 a/5 + 4 a/5 *) 7]/ 16 * 2 



Here 7 = a/2 [8 * 2 - 4 (1 - a/5) * + 15 - 5 a/5] • 
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TABLE VI. Moments (sum over paths) Mg through M30 
for Cqq with t = 1. 



Order(/) 


Moment (Mi) 


9 


—306 


10 


4274 


11 


—3080 


12 


31227 


13 


—29718 


14 


234559 


15 


—279100 


16 


1803375 


17 


—2572542 


18 


14149891 


19 


—23398880 


20 


113056535 


21 


—210843318 


22 


918114387 


23 


— 1887655172 


24 


7564926707 


25 


-16828070362 


26 


63140353799 


27 


-149626028160 


28 


532999985631 


29 


-1328522904154 


30 


4543918293899 



TABLE VII. Parameters a n and b n calculated from the moments Mi with t = 1 in the moment method for Ceo- 
in denominator are relatively prime to those in numerator and all numbers in square root are primes. 



Numbers 



n 

7T~ 
1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 

12 
13 
14 





_ l. 

3 
11 

69 

6633 
10925 
110900') 
2724600 
52107413 
(in42in7i;o 
66333080:517 



52324285077614 
603116478351886109 
1730319460378457102 
21011937073656177526717 
5<>J 7 I >2!.'«7i. ( ) I U..~,:i<)00009 
■ ■ ■ . i -.. i i - 

156471266768327126141170869441 
130318231516558926494126646931921 
26315692870966750872559328607407057 
284524265315564231563836760572796495 

" (i! )n -2 2 33 140426800959 1398461 790 72 3-~> ^3 1 
290400229158644153 
220686749524952341 




^ 72-23-239 
g^g V5- 19-23857 
rV23- 239 -41357 



2853633 
23 



119: 



31455736123 
7 



3- 5 -11 -31 -80657 -23857 
V3- 7- 17- 41357 -11492779 
~r V3 • 11 • 31 - 137 - 80657 - 15257198651 



>| ., r>() , (i()g3 V2 ■ 17 ■ 23 ■ 41 ■ 67 ■ 11492779 ■ 432553690463 
83175749 3 1391302 7 V137- 15257198651 ■ 3163866047896700891 
RRn f^ QR7nn „ 01 V7 ■ 13 ■ 41 ■ 67 ■ 4544717 ■ 57739553 ■ 432553690463 
5 • 13 • 17 • 89 • 3163866047896700891 



22()(iS(i74<)524<)52341 



1G 



TABLE VIII. Moments (sum over paths) A4g through A4i6 for Ceo with an arbitrary t. 



Order(7) Moment (Mi) 

9 -54t 4 - 18 P - 162 t 2 - 72 

10 £ 10 + 50 1 8 + 40 t 7 + 500 £ 6 + 300 1 5 + 1470 + 380 i 3 + 1280 t 2 + 254 

11 - 154 1 6 - 132 - 946 t 4 - 286 t 3 - 1232 i 2 - 330 

12 i 12 + 72 1 10 + 80 i 9 + 1 1 10 + 1080 t 7 + 5640 i 6 + 3384 + 10200 t 4 + 2520 i 3 + 6192 t 2 + 948 

13 -364 t s - 546 t 7 - 3848 i 6 - 2886 1 5 - 10088 i 4 - 2782 i 3 - 7774 1 2 - 1430 



i 14 + 98 1 12 + 140 i 11 + 2156 t 10 + 2940 t 9 + 16870 1 8 + 16464 t 7 + 
52444 i 6 + 30128 1 5 + 65016 i 4 + 15232 t 3 + 29456 1 2 + 3614 

-756 1 10 - 1680 t 9 - 12390 1 8 - 15750 1 7 - 55330 1 6 
-36936 i 5 - 84690 t 4 - 21430 1 3 - 44130 1 2 - 6008 

i 16 + 128 1 14 + 224 1 13 + 3808 1 12 + 6720 i 11 + 42560 1 10 + 58240 i 9 + 203576 i 8 + 
189696 1 7 + 431264 1 6 + 234464 t 5 + 392720 1 4 + 87168 i 3 + 138816 i 2 + 13990 



TABLE IX. Moments (sum over paths) for states \I + > 
and \I- > with t = 1. 



Order(l) < 7+|H'|J+ > < J_|H'|J_ > 

9 -312 -300 

10 4319 4231 

11 -3278 -2882 

12 32339 30115 

13 -33436 -26000 

14 252339 216779 

15 -225080 

16 1571823 



TABLE X. Parameters a n 
ments < 7±|H'|I± > with t ■ 



and b n computed from the mo- 
= 1. 





From < J + 


|H'|/+ > 


From < I- 


|H ! |/_ > 


n 


a n 




a„ 

















1 













21850 
1964421 

" 36832450 
40798 




^V3- 23 -137 -283 

855 



,-V2-3 - 13- 19 



_ i 

3 
11 
69 
1701 
4370 
11563 
44650 
342588 
22258025 
3661 



V2_ 

23 



|V2_E9_ 
5- 23^ 47 

Ti33-Vl7-47-89 
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TABLE XI. Moments (sum over paths) < I±\H l \I± > (I = 9, • • • , 16) for an arbitrary t. 



Order(Z) 



Moment (< J 4 



:h'|j+ » 



-54t* - 24?"- 162^ - 72 
t 1 " + 50 t 8 + 40 t 7 + 500 1 6 + 304 f° + 1490 1 4 + 400 i 3 + 1280 t 2 + 254 
- 154 1° - 176 1 5 - 990 1 4 - 396 1 3 - 1232 t 2 - 330 
t 12 + 72 i 10 + 80 t 9 + 1 1 10 t 8 + 1 104 t 7 + 5784 i 6 + 3648 + 10560 t 4 + 2840 i 3 + 6192 t 2 + 948 
-364 1 8 - 728 t 7 - 4212 i 6 - 4056 1 5 - 10920 i 4 - 3952 t 3 - 7774 f 2 - 1430 



9 

10 
11 
12 
13 



14 



i 14 + 98 f 12 + 140 f 11 + 2156 1 10 + 3024 t 9 + 17458 1 8 + 18032 1 7 + 
55944 i 6 + 35000 1 5 + 69020 i 4 + 18396 t 3 + 29456 1 2 + 3614 



15 



-756 t 10 - 2240 t 9 - 14070 1 8 - 22380 i 7 - 64470 t e 



-53736 i 5 - 94230 t 4 - 31100 1 3 - 44130 i 2 



6008 
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Order(Z) 



t 16 + 128 t 14 + 224 1 13 



3808 1 12 + 6944 i 11 



- 44352 r 



• 64512 f + 222168 i s + 



225056 1 7 + 481696 r + 292960 1° + 428144 t 4 + 112128 1 3 + 138816 1 2 + 13990 



Moment (< J_|H'|/_ >) 



-54 1* - 12 F- 162^ - 72 
£ 10 + 50 + 40 i 7 + 500 + 296 1 5 + 1450 1 4 + 360 1 3 + 1280 1 2 + 254 
-154 1 6 - 88 1 5 - 902 i 4 - 176 i 3 - 1232 1 2 - 330 
i 12 + 72 t w + 80 i 9 + 1110 1 8 + 1056 1 7 + 5596 1 6 + 3120 t h + 9840 i 4 + 2200 1 3 + 6192 1 2 + 948 
-364 - 364 1 7 - 3484 - 1716 1 5 - 9256 i 4 - 1612 1 3 - 7774 i 2 - 1430 



9 

10 
11 
12 
13 



14 



f 14 + 98 1 12 + 140 1 11 + 2156 f 1(1 + 2856 1 9 + 116282 f 8 + 14896 1 7 + 
48944 + 25256 t s + 61012 f 4 + 12068 i 3 + 29456 1 2 + 3614 



15 



-756 1 10 - 1120 i 9 - 10710 1 8 - 9120 1 7 - 46190 1 6 
-20136 £ 5 - 75150 1 4 - 11760 i 3 - 44130 1 2 - 6008 



16 



t 16 + 128 1 14 + 224 1 13 +3808 1 12 + 6496 1 11 + 40768 1 10 + 51968 1 9 + 184984 1 8 + 
154336 i 7 + 380832 t G + 175968 £ s + 357296 1 4 + 62208 1 3 + 138816 1 2 + 13990 



TABLE XII. Moments (sum over paths) < O |H'|0 O > and < (p \H l \ipo > for / = 1, • • • , 8. 



Order(Z) Moment (< O |H'|0 O >) 

1 ~+2 

2 t 2 + 4 

3 -4t 2 -8 

4 + 14t 2 + 16 

5 -6t 4 - 2i 3 - 40 i 2 - 32 

6 £ 6 + 30t 4 + 8i 3 + 110i 2 +64 

7 -8t fi - 6t 5 - 112 1 4 - (36 + 2T)t 3 - 282 1 2 - 128 

8 t 8 + 52 1 6 + 32 1 5 + (396 + AT) t 4 + (1 16 + 12P) i 3 + 708 1 2 + 256 

Order(l) Moment (< y |H'|y >) 

1 r+i 

2 t 2 + 2 1 + 2 

3 t 3 + 3 i 2 + 5 i + 3 

4 t 4 + 4 £ 3 + 10 1 2 + 10 1 + 5 

5 t 5 + 5t 4 + 16t 3 + 25t 2 + 20t + 8 

6 t 6 + 6t 5 + 24t 4 + 48t 3 + 60t 2 + 384 + 13 

7 t 7 + 7t 6 + 33 1 5 + 84 1 4 + (133 + P) t 3 + 133 1 2 + 71 1 + 21 

8 t 8 + 8 1 7 + 44 1 6 + 132 1 5 + (266 + 2T) t 4 + (336 + AT) t 3 + 284 1 2 + 130 1 + 34 
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FIG. 1. Flattened Ceo molecule obtained by stretching (a) a pentagon and (b) a hexagon, with the site labels used in cases 
A and B in the recursion method. 



FIG. 2. Site labels used in the path-integral method. Independent vertices are those with j running from 1 to 24. 

FIG. 3. The local density of states around ring-clusters : (a) a pentagon, (b) a hexagon. They are obtained by using the 
recursion method. Notice that around a pentagon (hexagon) ring the LDOS is large at low (high) energies. This is related to 
the fact that a pentagon (hexagon) has zero (three) double bonds and five (three) single bonds. 

FIG. 4. The local density of states around two opposite (i.e., antipodes) pentagons with (a) even parity, (b) odd parity, and 
two two opposite hexagons with (c) even parity, (d) odd parity. They are obtained by using the recursion method. Notice that 
around two antipode pentagon (hexagon) rings the LDOS is large at low (high) energies. 

FIG. 5. The local density of states on two opposite carbon atoms (antipodes) with (a) even parity, (b) odd parity. They are 
obtained by using the moment method. Notice that there are 7 even parity solutions in (a) and 8 different odd parity solutions 
in (b); for a total of 15 different energy levels (for t = 1). 

FIG. 6. The local density of states for a carbon atom obtained from the moment method with the respective degree of 
degeneracy shown above each peak. 
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